Asymmetric diffusion at the interfaces in multilayers 
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Nanoscale diffusion at the interfaces in multilayers plays a vital role in controlling their physical 
properties for a variety of applications. In the present work depth-dependent interdiffusion in 
a Si/Fe/Si trilayer has been studied with sub-nanometer depth resolution, using x ray standing 
waves. High depth-selectivity of the present technique allows one to measure diffusion at the two 
interfaces of Fe namely, Fe-on-Si and Si-on-Fe, independently, yielding an intriguing result that Fc 
diffusivity at the two interfaces is not symmetric. It is faster at the Fe-on-Si interface. While the 
values of activation energy at the two interfaces are comparable, the main difference is found in 
the pre-exponent factor suggesting different mechanisms of diffusion at the two interfaces. This 
apparently counter-intuitive result has been understood in terms of an asymmetric structure of the 
interfaces as revealed by depth selective conversion electron Mossbauer spectroscopy. A difference 
in the surface free energies of Fe and Si can lead to such differences in the structure of the two 
interfaces. 



Atomic diffusion is fundamental to many processes 
in material science such as microstructure development, 
non-martensitic phase transformation, stress relaxation 
etc. In multilayers, atomic diffusion at the interfaces 
plays a vital role in controlling their physical properties 
for a wide variety of applications. In x ray and neutron 
mirrors intentionally diffused interfaces have been used 
to reduce the higher order contamination In tunnel 
magnetoresistance multilayers, thermal annealing can in- 
crease magnetoresistance by orders of magnitude 2] . In 
giant magnetoresistance (GMR) multilayers, interdiffu- 
sion can significantly affect the GMR [3J. In spin valves 
with Mn based antiferromagnetic layer, Mn diffusion can 
seriously degrade the performance Co diffusion in 
Sm-Co/Fe exchange-spring magnet films is known to im- 
prove the exchange coupling [5J. While atomic diffusion 
in bulk solids is a widely studied and fairly well under- 
stood phenomenon, a reasonable understanding of the 
interfacial diffusion in multilayers has yet to emerge. Sev- 
eral factors like a steep concentration gradient at the 
interfaces, interfacial stresses and disorder may signifi- 
cantly modify the diffusion in multilayers. This has re- 
sulted in unexpected interfacial phenomenon like a non- 
parabolic shift of phase boundaries in the presence of 
strong composition dependence of diffusivity Q . 

X ray standing waves generated by total external re- 
flection of x rays from buffer layer Q, or in a multi- 
layer H 0] have been used for concentration profiling 
of various elements. In the present work, we exploit 
the depth selectivity of x ray standing waves for study- 
ing the depth dependent interdiffusion in a Si/Fe/Si tri- 
layer with sub-nanometer depth resolution. High depth 
sensitivity of the technique allows one to measure dif- 
fusion at the two interfaces of Fe namely Fe-on-Si and 
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Si-on-Fe independently. An intriguing finding of the 
present work is that the Fe diffusivity at the two in- 
terfaces is not symmetric; diffusion is faster at the Fe- 
on-Si interface. In order to understand this apparently 
counter- intuitive result, conversion electron Mossbauer 
spectroscopy (CEMS) has been used to study the inter- 
facial structure. It is found that there is a significant 
difference in the structure of the two interfaces, resulting 
in different diffusivities at the two interfaces. A differ- 
ence in the surface free energies of Fe and Si can lead to 
such difference in the structure of the two interfaces. 

The structure of the multilayer used for diffusion mea- 
surements is [W (2.0nm)/Si (3.1nm)]in / Si (3.8 nm)/ Fe 
(2.7nm)/Si (7.0 nm) (referred as SW_ML). The deposi- 
tion was done using ion beam sputtering in a vacuum 
chamber with a base pressure of lxlCP 7 mbar A 
broad beam Kaufman type ion source was used with Ar 
ions of energy IkeV and a beam current of 20 mA. The 
bottom [W (2.0nm)/Si (3.1nm)Jip multilayer is used to 
generate x ray standing waves i|. On the top of this 
multilayer Si (3.8nm)/Fe (2.7nm)/Si (7.0 nm) structure 
was deposited without breaking vacuum. A thickness of 
3.8 nm for the first Si layer is chosen in such a way that 
the Fe layer lies roughly midway between two antinodes 
of the x ray standing waves generated in W/Si multi- 
layer at the Bragg peak. This point is clear from the 
inset of Fig. [T] which shows the contour plot of x ray in- 
tensity as a function of depth and the scattering vector 
q. At q — 1.35 nm -1 which corresponds to the center of 
the Bragg peak of W/Si multilayer, Fe layer is midway 
between the two antinodes. 

Simultaneous x ray reflectivity (XRR) and x ray flu- 
orescence (XRF) measurements were done using Brukcr 
D8 diffractometer fitted with a Gobble mirror on the in- 
cident beam side in order to obtain a parallel monochro- 
matic beam of Cu Kcv radiation. Fluorescence spectrum 
was measured using a Kctek detector with an energy res- 
olution of 200 eV. Thermal annealing of multilayer struc- 
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FIG. 1: Representative XRR of SW-ML after various stages 
of annealing. The continuous curves represent the best fit 
to the experimental data. The numbers shown against each 
curve represent annealing temperature and the average inter- 
face roughness of W/Si multilayer respectively. Inset shows 
the contour plot of x ray intensity as a function of q and depth 
from the surface of the multilayer. The rectangle represents 
the position of Fe layer. 



ture was done in a vacuum of 10~ 6 mbar in order to 
induce interdiffusion at the interfaces of Fe and Si lay- 
ers. A separate study on the annealing behavior of W/Si 
multilayer showed that at least up to 623 K, the W/Si 
multilayer structure is stable with almost no change in 
its x ray reflectivity. Thus the thermal annealing in the 
present multilayer essentially induces interdiffusion at the 
interfaces of Fe and Si layers only. 

In order to study the Fe-on-Si and Si-on-Fe interfaces 
using CEMS, the following two multilayer structures 
were also prepared using ion beam sputtering: (i) sub- 
strate/ [Si(6nm)/ 57 Fe(2nm)/Fe(2nm)]i , and (ii) sub- 
strate/ [Si(6nm)/Fe(2nm)/ 57 Fe(2 nm)]io, designed as 
ML1 and ML2 respectively. In the first multilayer the 
2nm thick Fe layer lies on Fe-on-Si interface while in 
the second multilayer it lies on Si-on-Fe interface. Since 
Mossbauer measurements are sensitive only to 57 Fe, the 
first multilayer gives information preferentially about the 
Fe-on-Si interface, while the second one gives information 
about Si-on-Fe interface. The CEMS measurements were 
done using a Wissel Mossbauer spectrometer and a gas 
flow proportional counter with 95% He + 5% CH 4 . 

SW_ML was isochronally annealed at 473 K, 498 K, 
523K, 598K and 623K for lh each. Figure [1] gives some 
representative XRR of SW-ML after various stages of 
annealing. The XRR of this multilayer is dominated by 
that of W/Si multilayer mirror, with only small modula- 
tion in the q region below the Bragg peak attributable to 
the top Si/Fe/Si trilayer. Figure [2] shows Fe- fluorescence 
data of SW_ML measured simultaneously with XRR. 
The region around the Bragg peak in which x ray stand- 
ing waves are generated, is highlighted in the figure by 
shaded area. One may note that in the as-deposited sam- 
ple Fe-fluorescence exhibits two well defined peaks in the 
region where x ray standing waves are generated. The 
origin of these two peaks can be understood from the 
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FIG. 2: Fe-fluorescence of SW-ML after various stages of 
annealing. The continuous curves represent the best fit to 
the experimental data.The shaded area represents the region 
around Bragg peak where x ray standing waves are formed. 
The dashed curve represents the best fit to the experimental 
data of 498K annealed sample with the roughnesses of the 
two interfaces of Fe taken to be equal. 



contour plot of x ray intensity (inset of Fig [T]) . At the 
center of the Bragg peak, Fe layer lies roughly midway be- 
tween the two antinodes. However, as one moves away on 
either side of the Bragg peak, the antinodes get shifted 
resulting in partial overlapping of one of the antinodes 
with Fe layer, giving rise to a peak in Fe-fluorescence. 
The peak around q = 1.27 nm -1 corresponds to a situ- 
ation where one of the antinodes partially overlaps with 
Fe-on-Si interface while the peak at q = 1.41 nm -1 occurs 
as a result of partial overlap of an antinode at Si-on-Fe 
interface. With thermal annealing both these peaks get 
broadened and their intensities get modified. However, 
changes occuring in the two peaks are quite different, sug- 
gesting that the two interfaces get modified differently 
with thermal annealing. Simultaneous fitting of XRR 
and Fe-fluorescence data has been done using Parratt's 
formalism The average roughness of the interfaces 
of W/Si multilayer owjSi-, an d the roughnesses of the two 
interfaces of Fe, UFe-on-Si an d crsi-on-Fe were the only 
parameters which were varied as a function of annealing 
temperature. While the XRR data is mainly sensitive to 
the changes in <7w/Sii the Fe-fluorescence is affected by 
the changes in the structure of the interfaces of Fe layer. 
One finds that with thermal annealing up to 623 K, there 
are only minor changes in the W/Si multilayer with the 
interface roughness going from 0.6 nm to 0.85 nm. The 
roughnesses of the Fe-on-Si and Si-on-Fe interfaces as ob- 
tained from the fitting of the fluorescence and XRR data 
are shown in Fig. O One may note that the roughness of 
Fe-on-Si interface increases at a much faster rate as com- 
pared to that of Si-on-Fe interface. It may be mentioned 
that for the thermal annealing upto 473 K, the width of 
Fe concentration profile is much smaller than the separa- 
tion between two antinodes. Therefore, from the fluores- 
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FIG. 3: Arrhenius plot of Fe diffusivity at the two interfaces. 



cence data it is difficult to estimate the roughnesses of the 
two interfaces separately. However, at 498K and above 
roughnesses of the two interfaces can be determined in- 
dividually with good reliability. In order to demonstrate 
this point, the fitting of 498 K data obtained by taking 
the roughnesses of the two interfaces to be equal is also 
shown in Fig. [2] (dashed curve). The best fitted curve 
clearly deviates from the experimental data. 

The variation in the roughnesses of the two interfaces 
with the thermal annealing has been used to estimate the 
diffusivity of Fe at the two interfaces using the relation: 
D{T 2 ) = [cr 2 (T 2 ) - cr 2 (Ti)]/2£,where D(T 2 ) is the diffu- 
sion coefficient at temperature T 2 , cr(Tx) and u(T 2 ) are 
the roughnesses of an interface before and after annealing 
at temperature T 2 for time t. Figure. [3] gives the Arrhe- 
nius plot of InD versus l/T at the two interfaces of Fe. 
One finds that there is a significant difference in Fe dif- 
fusivities at the two interfaces. A linear fit to the experi- 
mental data as expected from the Arrhenius temperature 
dependence of the diffusivity D(T) = Doexp(—E / 'KbT) : 
yields the activation energy E as well as the pre-exponent 
Dq for the diffusion. The experimentally obtained values 
are shown in Fig. [3] One finds that the activation en- 
ergies for diffusion at the two interfaces are comparable 
within experimental errors. However, the pre-exponent 
Dq is significantly high for Fe-on-Si interface as compared 
to that for Si-on-Fe interface, resulting in a significantly 
high diffusivity at Fe-on-Si interface in the temperature 
range studied in the present experiment. 

At first glance the above results appear to be counter 
intuitive as both interfaces ought to be identical, having 
Fe layer on one side and Si layers on the other. In order 
to understand the possible reason for this difference in 
the diffusivity at the two interfaces, a detailed study of 
the structure of the two interfaces has been done using 
CEMS on the samples ML1 and ML2. As discussed in 
the experimental section CEMS of ML1 is sensitive to Fe- 
on-Si interface while that of ML2 is sensitive to Si-on-Fe 
interface. 

Figure[4ja) shows Mossbauer spectra of ML1 and ML2 
in as-deposited state. The spectra of both the specimens 
are fitted with three overlapping components: (1) a sharp 
sextet with hyperfine field about 33 T, (2) a broad mag- 
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FIG. 4: (a) Conversion electron Mossbauer spectra of as- 
deposited ML1 and ML2,(b) The corresponding hyperfine 
field distribution of the broad magnetic sextet. 



TABLE I: Results of fitting of Mossbauer spectra of ML1 
and ML2 as described in the text. S, A and B^f represent 
isomer shift, quadrupole splitting and average hyperfine field 
respectively. Isomer shift 5 is with respect to a-Fe. 



Sample 


ML1 


ML 2 


Sharp Sextet B h/ (T) 
Area (%) 


31.34±0.08 

28±6 


32.32±0.06 
32±5 


Broad Sextet B h/ (T) 
Area (%) 


26.11±0.32 

46±6 


25.74±0.21 
50±3 


S (mm/s) 
Doublet A (mm/s) 
Area (%) 


0.16±0.01 
0.68±0.02 
26±3 


0.24±0.01 
0.65±0.03 
18±3 



netic component, having a distribution of hyperfine fields 
and (3) a non-magnetic doublet. The sharp sextet rep- 
resents the bulk of a-Fe, while the broad magnetic com- 
ponent and the doublet represent the Fe atoms in the 
interfacial region. Further it is known that in Fe-Si alloy, 
if iron concentration is less than 50%, it becomes non- 
magnetic [T3|. Therefore, the area under the doublet rep- 
resents the fraction of Fe atoms in the interfacial region 
having iron concentration less than or equal to 50%. Re- 
sults of fitting are given in table (1). One may note that 
the area under the sharp sextet in both the specimens is 
about 30%, within experimental errors. However, the rel- 
ative areas of broad magnetic component and the doublet 
as well as the shape of the hyperfine field distribution of 
broad magnetic component [Fig. b)] are very different 
in the two specimens. This suggests that the structure 
of the two interfaces namely Si-on-Fe and Fe-on-Si inter- 
face may be different. A rough estimate of the interface 
roughness as obtained from the width of the intermixed 
region comes out to be 0.78 nm, which is in agreement 
with that obtained from x ray measurements. From table 
(1), one finds that at Si-on-Fe interface, the fraction of 
Fe atoms in the doublet is 18%, which is about 23% of 
the total Fe atoms in the intermixed region. This area 
fraction agrees very well with that expected for an error 
function concentration profile [Fig. [S]- Thus, Mossbauer 
measurements suggest that the concentration profile at 
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FIG. 5: Schematic diagram of concentration profile of Fe 
layer. The regions corresponding to sharp sextet, broad sex- 
tet and doublet are shown in different shades. The dashed 
rectangle shows the position of ideal interfaces. 

Si-on-Fe interface is an error function. 

From Mossbauer spectrum of specimen ML1 one finds 
that in this spectrum, the area of doublet is 27%, which 
is significantly more than what one expects for an er- 
ror function concentration profile. Further a compari- 
son of hyperfine field distributions at the two interfaces 
shows that at Si-on-Fe interface, there is a continuous 
distribution of hyperfine field values starting from 32 T 
down to 10 T while at Fe-on-Si interface, low field com- 
ponents are missing. The above two differences can be 
understood, if one assumes the concentration profile at 
Fe-on-Si interface to be as shown schematically in Fig. [5j 
where the area from the low field magnetic component 
gets transferred to the non-magnetic component. This 
suggests that at Fe-on-Si interface, there is an interlayer 
of Fei_ x Sia;, resulting in a plateau in a concentration 
profile. The isomer shift and quadrupole values of the 
doublet of ML1 match very well with those of FeSi2, sug- 
gesting that the composition of interlayer is FeSi2 (l3j |. 
Thus, at Fe-on-Si interface diffusion occurs via the FeSi2 
interlayer, while at Si-on-Fe interface diffusion occurs via 
bcc Fe (Si) phase. In literature the diffusivity data for 
self-diffusion in stiochometric DO3 phase as well as inter- 
metallic compounds of FeSi are given. However, it is not 
meaningful to compare the diffusivities obtained in the 
present work at the two interfaces with bulk diffusivities, 



since it is known that the multilayer diffusivities at the 
interfaces are very different because of high concentration 
of defects and possible concentration gradient [3]. 

The difference in the structure of the two interfaces 
can be understood in terms of the difference in the sur- 
face free energies of Fe(2.9Jm~ 2 ) and Si (1.2 Jm~ 2 ) [l5[. 
During the deposition of Fe on Si, the surface free energy 
of Si being lower, Si atoms try to move to the surface 
guided by the chemical driving force. This would lead to 
a stronger mixing at the interface and a possible forma- 
tion of FeSi2 compound. On the other hand, during the 
deposition of Si on Fe, no such chemical driving force ex- 
ists, therefore the intermixing at Si-on-Fe interface would 
take place as a result of random thermal motions only 
and hence concentration profile is expected to be an er- 
ror function. 

In conclusion, x ray standing wave technique has been 
used to get concentration profile of Fe layer in Si/Fe/Si 
trilayer. The precision of this technique is sufficient to 
differentiate between the two interfaces of the Fe layer. 
This allows one to study interdiffusion at the two inter- 
faces namely Fe-on-Si and Si-on-Fe independently. Inter- 
estingly the diffusivities at the two interfaces are signif- 
icantly different. This seemingly counter- intuitive result 
can be understood in terms of a difference in the struc- 
ture of the two interfaces in the as-deposited film itself. 
CEMS measurements show that while at Si-on-Fe inter- 
face, the Fe concentration profile is an error function, at 
Fe-on-Si interface an interlayer exists with approximate 
composition of FeSi2- This difference in the structure 
of the two interfaces is the cause of different diffusivi- 
ties at the two interfaces. Besides being of fundamental 
importance in understanding the interfacial diffusion at 
nanometer scale, present results have important impli- 
cations on the use of controlled thermal annealing for 
tailoring the properties of multilayers for a wide variety 
of applications. 
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